RU-NHETC-2012-12 
UTTG-09-12 



Di-Jet Extinction from Non-Perturbative 
Quantum Gravity Effects 



Can Kilic 1 , Amitabh Lath 2 , Keith Rose 2 , Scott Thomas 2 

1 Theory Group, Department of Physics and Texas Cosmology Center 
The University of Texas at Austin 
Austin, TX 78712 

2 Department of Physics 

Rutgers University 
Piscataway, NJ 08854 



Abstract 

We study a novel signature of TeV scale quantum gravity that manifests itself as an 
extinction of hard short distance scattering in QCD processes. The extinction behav- 
ior is due to the predominance of high-entropy intermediate states of the underlying 
quantum gravity theory. We model extinction using a large damping Veneziano form- 
factor modification of QCD scattering amplitudes that suppresses high pt scattering. 
We propose and demonstrate the potential of an LHC search for extinction, with a 
possible reach for the string scale as high as 3 TeV with 7 TeV LHC collision data, 
and up to 5 TeV from high-statistics 13 TeV data. 



1 Introduction 



The LHC, currently operating at 8 TeV center of mass energy, and planned to operate at 
13 or 14 TeV after the shutdown period beginning at the end of 2012, is exploring new 
frontiers in short distance physics. Among the most dramatic outcomes of probing physics 
beyond the TeV scale is the possibility of large extra dimensions [1] where the fundamental 
Planck scale is exponentially reduced, which may make it possible for the LHC to access 
states associated with the fundamental properties of gravity at the shortest scales. A 
reduced Planck scale in the TeV regime can also be realized in theories with warped extra 
dimensions [2]. Various potential signatures of such a low Planck scale have been considered 
in several guises [31 HJ El E] , including in the form of string Regge resonances [3 El [9] when 
the underlying string theory is weakly coupled, as well as the formation of microscopic 
black holes JTU1 [HI H21 EE31 EH EES] in the strongly coupled regime. In this work we wish to 
emphasize a very different type of signature of strongly coupled TeV scale quantum gravity, 
which may give the earliest indication of the exotic underlying theory. 

Specifically, in realizations of TeV scale gravity based on string theory in a strong 
coupling regime, a novel feature at center of mass energies comparable to the string scale 
is the contribution to scattering processes from the exchange of high entropy states. This 
feature is not specific to any single underlying string theory model but is rather robust. 
Due to the high entropy of these intermediate states, processes with a small number of 
final state particles become suppressed [HI [T7]. As we will show later in this paper, this 
effect manifests itself immediately at the string scale, while the production of black holes is 
expected to be realized only at significantly higher energies. Thus the earliest observable in 
such a scenario may be the "extinction" of all 2 — > 2 processes with center of mass energies 
of roughly the string scale. In the context of the LHC, this will be observable as an 
exponential drop in the high-pr QCD di-jet production. It should also be emphasized that 
in collider studies of black hole production, the models that are used exhibit a discontinuity 
as the cross section for black hole formation is taken to jump abruptly from zero to a large 
geometric cross section at the threshold energy. Extinction on the other hand can be 
modeled in a way such that the turn-on of the suppression is gradual, and the transition 
from QCD-like behavior to the extinction regime is smooth. 

In the following study, we will model the effects arising from the underlying string theory 
description as Veneziano type form-factors that modify all 2 — > 2 scattering amplitudes in a 
universal manner. Based on this model, we will then derive the form-factor modified QCD 
scattering amplitudes for various initial and final states. Apart from studying these mod- 
ifications analytically, we also implement the modified amplitudes into a matrix-element 
based Monte Carlo event generation tool and thus perform a realistic collider study of the 
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extinction phenomenon, and assess the sensitivity of the LHC to a string scale in the TeV 
regime with early data as well as with large integrated luminosity after the energy upgrade. 



2 Modeling Extinction 

2.1 Damped Veneziano Scattering Amplitudes 

The goal of this section is to model the effects due to the exchange of high-entropy states 
by a phenomenological low energy description, and then to study the most salient features 
quantitatively. The idea behind the model is to map the color-ordered helicity amplitudes 
for any 2 — > 2 QCD process on to the corresponding string theory amplitude. 
Each string theory amplitude includes the Veneziano form factor 

j r / v r(i -x)r(i -y) 

v <™ )= r ( i -*-„)) (1) 

where x, y are kinematic invariants derived from the momenta associated with neighboring 
operators (as defined below). 

To match this behavior, the color ordered helicity amplitudes for external states with 
momenta Pi,P2,P3,Pa are taken to be those of QCD modifed by the Veneziano form factor 
in adjacent kinematic channels in the color ordered amplitude 

A(pi,P2,P3,Pa) = A(pi,P2,Ps,Pa)qpd ^(^12,^23) (2) 
where xab = (pa+Pb) 2 = ^Pa'Pb with all momenta incoming and 

V {x ,y ) = V^^,^±^). (3) 

Here M is the mass scale for the form factor modification, and a is a dimensionless damping 
parameter discussed below. The form factor approaches unity for small invariant momenta, 
V(0, 0) = 1, so the scattering amplitudes smoothly reduce to QCD at fixed momenta in 
the limit of large mass scale 

lim A(pi, p 2 ,Ps, Pa) = A(p 1 ,p 2 ,p 3 , Pa) qcb (4) 

M— >oo 

Let us now derive the form-factor modified amplitudes for all QCD processes, organized 
by different choices of external states. 
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2.1.1 Two Quarks of One Flavor and Two Quarks of Another Flavor 

The unique color ordered helicity amplitudes with two quarks of one flavor and two quarks 
of another flavor may be written 

AQi\ 4\ Qs ±k , 4?) = A (Qi\ ti J , (T a )/(T a ) e k (5) 

with others related by charge conjugation or parity, where throughout helicity and momen- 
tum indicated in helicity amplitudes are all incoming. The . . . are (anti)fundamental 
indices, a,b, . . . are adjoint indices, and T a are fundamental generators of SU(3)c with 
normalization Tr(T a T b ) = 5 ab . The QCD helicity amplitudes with color factors stripped off, 
written in terms of twistor products are related by permutations to a single amplitude 

A(n-* n+i n'- k n'+ l \ - n 2 ^ ^ (fC\ 

A K<ll , ?2 > % > ?4 JQCD - 9 {12) [21] ^ ' 

At?!-* n+j a'+ k n'~ e \ - n 2 ^ j] ^) / ? n 

a {Qi ,9a >?3 ,94 Jqcd - 9 {12) [21] ' ' 

where the twistor products are defined in terms of contractions of the wave function spinors 
by [AB] = x a (PA)Xa(PB) and (AB) = xI{pa)x*"{pb) and satisfy (AB) — —(BA), (AB)* = 
[BA\. These color ordered helicity amplitudes may be averaged (summed) over intial (final) 
state color and helicity configurations to give the averaged squared matrix elements for 
2 — > 2 scattering processes involving two quarks of one flavor and two of another flavor 
with the form factor modification 

(\A( qi q 2 q^)\ 2 ) = g 4 j^-j^ \V(s,t)\ 2 (8) 
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\V{t,u)V (9) 

(1-4(9^2 MDI 2 ) = g 4 '^ (s ' t)|2 (10) 

(\AM M 4 )| 2 > = 9 4 \ S -^- \V{t,u)\ 2 (11) 

where (AB)[BA]=xab, and s = Xi 2 = x 34 , t = Xi 3 = x 2 4, and m = xi 4 = x 2 3 are the Mandel- 
stam variables for the scattering order 12—7-34 indicated in the squared matrix elements. 
Each scattering matrix element squared flHl fTTj) includes a sum over two independent helic- 
ity channels (not related by parity) of single helicity amplitudes squared in a single color 
channel. 
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2.1.2 Four Quarks of the Same Flavor 



The color ordered helicity amplitudes with four quarks of the same flavor with either all the 
chiralities the same or two of one chirality and two of the other may be written respectively 

4\ %\ if) = a (qi\ qt\ %\ if) cn/(n/ (12) 

A(qr, Q2 J , qt\ qf) = df k , <l¥) (T a )/(T a ) e k (13) 

with others related by charge conjugation or parity. The QCD helicity amplitudes in this 
case with the color factor stripped off are identical to the amplitudes and ([7]). These 
color ordered helicity amplitudes may be averaged (summed) over intial (final) state color 
and helicity configurations to give the the averaged matrix elements squared for 2—7-2 
scattering processes involving four quarks of a single flavor with the form factor modification 

4 

(\A(qiq 2 ^<M 4 )| 2 ) =9 4 - 

4 

(\A( qi q 2 ^ q 3 q4)\ 2 ) = Q A - 

4 

{\A(M2 ^M±)\ 2 ) = 9 g 

Each scattering matrix element squared (I14H16P contains two types of terms. The first 
type for two quarks of one chirality and two of the other chirality includes a sum over two 
independent helicity channels (not related by parity) of single helicity amplitudes squared 
in a single color channel. The second type for four quarks of the same chirality includes a 
sum over the same two independent helicity channels of the square of a coherent sum of 
two helicity amplitudes in different color channels. 
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2.1.3 Four Gluons 

The unique four gluon color ordered helicity amplitudes may be written 

A9i a , 9t\ gf c , 9T) = Afa, g+, gf, gf) Tr(T a T b T c T d ) (17) 

A9i a , 92 b , 9t c , 9t d ) = Afa, g;,gt g+) T r (T a T b T c T d ) (18) 

with others related by parity. The QCD helicity amplitudes with color factors stripped off 
are all related by permutations to the Parke- Taylor amplitude [T8] 



(14) 4 

A( gi , gt, gf, <7 4 )qcd = ^ (12 )( 2 3>(34>(41> (19) 
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(13) 

A( ^'^ + '^'^ +)qcd= ^(12)(23)(34)(41) (20) 

(12) 4 

A(gi, g 2 ,gt, gthcv = g 2 (12) (23) (34) (41) (21) 

These color ordered helicity amplitudes may be averaged (summed) over intial (final) state 
color and helicity configurations to give the the averaged matrix elements squared for 2 — > 2 
scattering processes involving four gluons with the form factor modification 



[\A( gm gm ) I 2 ) = g 4 \ + - + ^) 



s 2 \V(t,u)\ 2 + t 2 \V(s,u)\ 2 + u 2 \V(s,t)\' 



4 , ~ ~ ,2 

- -=\s V(t,u)+t V(s,u)+u V(s,t)\ 



(22) 



27' 

The scattering matrix element squared (|2"2"|) includes a sum over three independent helicity 
channels (not related by parity) of the square of a coherent sum of three helicity amplitudes 
each in a different color channel. 

2.1.4 Two Quarks and two Gluons 

The unique color ordered helicity amplitudes for two quarks of a given flavor and chirality 
and two gluons may be written 

A(h\ 4\ gf\ gf) = Afo, g 2 + , gf, gf) (T a T b )/ (23) 

with others related by charge conjugation or parity The QCD helicity amplitudes with 
color factors stripped off are all related by permutations to a single amplitude 

+ + 2 (14) 3 (24) 

A[q x ,q 2 ,g 3 ,g 4 ) QCD - g (12)(2 3 )(34)(41) ( 24 ) 

(13) 3 (23) 

A(qT, q t,g;,gthc» = g 2 (12)(23)(34)(41) (25) 

These color ordered helicity amplitudes may be averaged (summed) over intial (final) state 
color and helicity configurations to give the the averaged matrix elements squared for 2 — > 2 
scattering processes involving two quarks of a given flavor and two gluons with the form 
factor modification 



(\A{giq 2 -> gzgi)\ 2 ) = g' 
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[\-A(gig 2 -»■ gs^) I 2 ) =g 
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^\V(s,t)\ 2 + -\V(s,u)\ 2 --Re(V(s,t)V*(s,u)) 
t u 4 v 

^\V(s,t)\ 2 + -\V(s, U )\ 2 --Re{V(s,t)V*(s,u)) 
t u 4 ' 



(26) 
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28) 



(29) 

Each scattering matrix element squared (I26H29P includes a sum over two independent he- 
licity channels (not related by parity) of the square of a coherent sum of two helicity 
amplitudes in different color channels. 

2.2 Quantitative Behavior of the Form Factor 

Unlike the case of the form factor with a small imaginary part, which leads to resonance- 
like behavior corresponding to string Regge excitations, in the case of complex argument, 
the form factor leads to exponential suppression. As we have discussed in the introduc- 
tion, this is the correct behavior we expect from the presence of high-entropy intermediate 
states. The asymptotic behavior of the form factor can be easily illustrated using Stirling's 
approximation for the gamma-function 



T(z) « V2ttz— e~ z . (30) 



Using the form-factor modified amplitudes listed in section 12.14 we illustrate the effect 
of extinction for several processes for central scattering (u — t — — |) in figure [TJ What is 
plotted here is the suppression relative to QCD, given by 

(31) 

\M? QCD (31) 

It is straightforward to check that for exactly central scattering, the suppression of all QCD 
processes listed in the appendix coincide with one of the four curves shown in figure [U It 
is easy to see on this figure that the suppression becomes very significant as soon as the 
center-of-mass energy of the collision becomes larger than the scale M. 

Of course, not all collisions at the LHC are central, so let us study the behavior of the 
form factor over the whole kinematic region. This will in fact point the way for choosing 
a good kinematic variable in which to bin the data and obtain good sensitivity to the 
extinction behavior. The 2 — ?- 2 process that dominates the cross section at high energy 
(with or without the form-factor modification) is qg — > qg so the behavior of the form factor 
for this process is a good approximation to the behavior for the sum of all 2 — > 2 processes. 
The suppression of qg — > qg as a function of scattering angle and center-of-mass energy 
is shown in figure |2j The black contours correspond to the log-base- 10 of the suppression 
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Figure 1: Effect of the Veneziano form factor extinction with a 
and gluon scattering probabilities for central scattering u — t- 
center-of-mass energy. 



1 on various 2—7-2 quark 
-s/2 as a function of the 



factor as defined in eq. ( |3T|) . as do the negative numbers labeling these contours, i.e. the 
contour labeled shows where the suppression is equal to 0.1 etc. We overlay on top of 
this the contours of constant p^, illustrated with the red contours, and the positive numbers 
labeling them in units of M 2 . We see that there is relatively good alignment between the 
contours of constant suppression and those of constant pt- This indicates that pt will be 
a good variable for binning the data, such that one is maximally sensitive to the effects of 
extinction. 

3 LHC Analysis 

3.1 Monte Carlo Implementation 

In order to study the effects of extinction at the LHC, we incorporate the form-factor- 
modified amplitudes into Pythia 6.4 |19j . Specifically, we modify the hard-scattering ma- 
trix elements occurring in the subroutine PYSGQC. With this modification we are able 
to generate events corresponding to the extinction scenario, and interface them to the 
(unmodified) parton shower and hadronization subroutines of Pythia. 

Since extinction is primarily an effect on large-angle scattering, we only generate hard 
scattering events within a pseudorapidity range of \r}\ < 1.5. Furthermore, in order to 
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Figure 2: Effect of the Veneziano form factor extinction with a = 1 on the qg — > qg 
scattering probability as a function of the center of mass scattering angle and center of 
mass energy. Black contours indicate \A(qg —> qg)\ 2 /\A(qg —> <7<7)|qcd labeled in log 10 
units. Red contours indicate p^/M 2 = (s/4)(l — cos 2 8)/M 2 in equal intervals. 

achieve high statistics in the large center-of-mass energy regime, we generate events in bins 
of Vs and combine the samples weighted by their relative cross sections. 

We also pass the resulting events through PGS [20j in order to account for detector 
effects such as acceptance and jet energy resolution. We use the CMS parameter set of 
PGS, with a jet definition based on a cone algorithm with AR = 0.7. 

3.2 Search Strategy 

Having demonstrated in section 12.21 that the form-factor suppression is well aligned with 
contours of p?, we will use histograms of pt in order to study the effects of extinction. 
Specifically, as our event selection criteria we demand that there are two hard jets in the 
event, both with px > 200 GeV and \rj\ < 1.5. 

Our analysis will be focused on evaluating the statistical significance with which models 
of extinction with various choices for the scale M can be distinguished from QCD. Clearly 
this should not be taken as a realistic estimate of the reach until various sources of sys- 
tematic uncertainty can also be included, such as PDF uncertainties and modifications of 
differential cross sections due to next to leading order effects. However, higher order cal- 
culations for 2 —7- 2 processes in the presence of quantum gravity effects are not currently 
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Figure 3: The differential cross section as a function of the p T of the hardest jet at 7 TeV 
after the selection cuts described in the text. 



available, and the same is also true for even leading order calculations of 2 — > 3 or higher 
multiplicity matrix elements. Therefore, currently an apples-to-apples comparison of QCD 
to extinction can only be performed at leading order. A comparison of extinction based 
on leading order matrix elements vs. next to leading order QCD is not recommended, as 
the extinction matrix elements are direct modifications of the QCD matrix elements by the 
form factor, therefore it would not be possible to model NLO QCD and extinction based on 
LO QCD such that they agree at low energies. In any case, since we are unable to account 
for such systematic uncertainties, our results should not be taken as precise statements 
about the LHC reach of an extinction search, but rather as a preliminary check that an 
extinction search is at least viable, even if the significance may be reduced by systematic 
effects. It should be noted that the log-likelihood analysis we will describe shortly relies 
on shape information only, therefore it is insensitive to at least certain types of systematic 
effects such as the scaling of the cross section by a constant k-factor. 

We thus proceed to assess the sensitivity of the LHC to extinction with the full 7 TeV 
dataset of 5 fb _1 , and with a future high-statistics sample of 100 fb _1 at 13 TeV. For the 
7 TeV study, we generate QCD events in bins of 300 GeV, starting at = 200 GeV, 
each bin containing 50000 events, which correspond to high statistics at the high end of 
the pt spectrum where we will perform our statistical analysis. We use the same V§ bins 
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Figure 4: The differential cross section as a function of the of the hardest jet at 13 TeV 
after the selection cuts described in the text. 



to simulate extinction with the scale M chosen at 1 TeV, 1.5 TeV, 2 TeV and 3 TeV, with 
25000 events in each bin, which is very high statistics considering the sharp decrease of 
extinction cross sections at high p?. The px spectrum of the hardest jet in each event are 
shown in figure [3j 

For the 13 TeV samples, we use bins of of 500 GeV starting at 500 GeV, with 100000 
events in each bin for QCD, and 25000 events per bin for extinction samples with a string 
scale of 2 TeV, 3 TeV, 4 TeV and 5 TeV. The pt spectrum of the hardest jet in each event 
is shown in figure HI 



3.3 Results 

We proceed to estimate the statistical significance with which the extinction scenario can 
be distinguished from QCD by setting up a log-likelihood analysis as follows. We divide 
the pt spectra shown in figures [3] and 0] in 100 GeV bins above pt = 500 GeV for the 7 TeV 
data and above p? = 1.5 TeV for the 13 TeV data. Using the cross section information in 
these bins we define the probability distributions for both the QCD hypothesis 

QCD 

P? CD = T^od (32) 

2-/ bins 
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as well as the extinction hypothesis defined in the same way. Note that the overall normal- 
ization cancels out such that the pi are sensitive to shape only. 

For each of the two competing hypotheses (QCD and extinction with a specific choice of 
the scale M) we perform 50000 pseudoexperiments based on these probability distributions, 
where the number of events per pseudoexperiment is chosen with a Poisson distribution 
around the mean number of expected events for that hypothesis. For each pseudoexperi- 
ment, the number of events in each bin is then used to define the log-likelihood ratio 



past the point where the two distributions intersect, we extract the p- value with which 
the two hypotheses can be distinguished. In both cases we find p-values that are below 
5.7 x 10~ 7 which corresponds to 5<r significance. Of course, as we mentioned before, in the 
absence of a more realistic treatment of systematic uncertainties, this result should only be 
taken as a first check for the plausibility of such a search. 



We have pointed to a novel signature of low scale quantum gravity, which manifests itself 
as a sharp decrease in all scattering processes with large transverse momenta. We have 
modeled this phenomenon by using form-factor modified QCD amplitudes based on un- 
derlying string theory considerations, where all scattering amplitudes are proportional to 
the Veneziano form factor. In order to model the behavior of the underlying theory in 
a regime of strong coupling, we have evaluated the form factors with large values of the 
damping parameter a, thus giving rise to the expected behavior of exponential suppression 
at large pr- We have performed a collider analysis of the extinction in order to assess 
the LHC sensitivity to such a scenario. Using a log-likelihood analysis based on statistical 
uncertainties only, we expect that a string scale of 3 TeV may be distinguished from QCD 
using the 7 TeV data sample, and with 100 fb _1 of integrated luminosity at 13 TeV this 
can increase to a string scale of 5 TeV. 

While future work is needed in order to include the effects of systematic uncertainties 
into this analysis and thus obtain a realistic estimate for the reach, this study shows that 
a search for extinction looks promising as a first indication of low scale quantum gravity at 




(33) 



Integrating the probability under the log-likelihood distributions of the two hypotheses 



4 Conclusions 



the LHC. 
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